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Abstract 


The  image  method  is  defined  and  all  regions,  equations,  and 
combinations  of  boundary  conditions  to  which  it  applies  are  de- 
termined. 
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I ,   Introduction 

The   ima^e   method  is   a  method  for  constructing  a  Green's   function  for 
a  part   of   space   bounded  by  planes    in  terms   of   the   corresponding  Green's 
function  for   the   full  space.      The  purpose   of   this   note   is    the  determination 
of  all   regions,   equations   and  boundary  conditions    to  v.'hich   this  method   is 
applicable. 

We   begin  by  considering  a  linear  second-order  differential  operator  L 
and  its   formal  adjoint  M,   both  defined   in  a  full   n-diniensional  Jiuclidean 
space.     We   may  define   a  Green's  function  of   the   operator  L  for  any  suffi- 
ciently regular  region  D  and  a  homogeneous    boundary  condition     on  the 
boundarj'  of  D,      This    is   a  function  G,r,r')    of   two  points   in  D  which  satisfies 
the   ho.T.ogenous   equation  MG  =  0  with  respect   to   Its   second  argument   r'   at  all 
points   in  D  except   r'   =  r,  at  whicn  point   it  has   an  appropriate  singularity. 
In  addition  G  satisfies    the  prescribed  boundary  condition.      If   the   region  D 
extends    to   infinity,   an  appropriate  homogenous   condition  must  also  be   imposed 
there. 

We  will  now  describe   the   imago   method  for  the   case    in  which  D   is   bounded 
by  hyperplanea   and   the   normal  derivative  vanishes   on  the   boundary  of  D, 
i'irst,  suppose  G   (r,r')   is   that  Green's  function  of  the  operator  L  for   the 
full   space  vrhich  satisfies   given  conditions   at  infinity.      Let  S(r)    be   the 
set   of   images   of  r   in   the   bounding  planes    of  D,    plus    the   images   of   these 
images    in  the  bounding  planes   of  D,    etc.     Now  we  construct   the   function 


(1)  0(r,r')  =  ^        G  *(r".  r'). 

r"£S(r) 


It   is  clear  that  S(r)    is  sjTiunetric    in  every  bounding  plane  of  C  and  therefore 
that  the  normal  derivative  of  G(r,r')  vanishes   on  the  boundaries  of  D,     For 
certain  operators   L,   each  term  of   the  sum  in   (l)    is   also  a  solution  of  M, 
Thus   for  these  operators,   G  will  be   the  Green's  function  for  D,   aissuming  the 
sum  converges,   and  provided   that   G  has   only  one   singularity  in  D.      This  will 
be    the  case   if   S(r)    has   only  one   point,   r   itself,    in  D, 

We  can  now  formulate  oxir  problems.      The  first  problem   is    to  find  those 
regions  D  with  the  property  that   the  set   S(r)    has   only  one  point,   r   itself, 
in  D,  for  every  point   r   in  D,      It   is   clear   that    these  are   the   only  regions 
for  which  the   image  method  described  above   is  applicable.      The   second  problem 
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is    that   of  finding,   for  each  such  region  D,   all   tnose   operators   L  for  vnich 
eacn  term  in   \1)    is   a  solution  of  MG  =  0 ,      It   is   also   clear   that  for  each 
adir.issible   region  E,    these   are    the   only  operators   for  which   tne   ims.ge   method 

just   described  applies.      /ne   third  problem  ccncernf    the   extension  of   the 

SG 
image   method  to  the  boundary'   conditions   G  =  0   and  r—  +  hG  -  0     and  coiiihina- 

"  on 

tions      of   ti-ese   conditions,      This   problem  will    be  f'orther   explained   in 
■iectiou  ili, 

II.  Adrcissible  Regions 

A  region  I)   for  wnich  the  image  method  is  applicable  we  will  call  admis- 
sible.  Let  E.  denote  a  reflection  in  side  i  of  an  aamissiole  region  D,   The 
1 

set  S(r)  of  images  of  ariy  point  r  in  D  is  the  set  of  points  obtained  by- 
applying  to  r  all  the  reflections  in  tne  group  generated  by  the  H  .   If  D  is 
admissible,  ii(r)  has  only  one  point,  r  itself,  in  C,  i3ut  this  implies  that 
tne  group  generated  by  the  H.  is  discrete  and  that  D  is  a  fundamental  doDiain 
of  tnis  group. 

However,  not  every  fundamental  domain  of  a  discrete  ^ro^P  of  reflections 
is  an  admissible  region  for  the  image  metnod.  In  fact,  the  angle  between  any 
two  planes  boxinding  the  fiandamental  domain  of  a  discrete  group  of  reflections 

is  of  tae  form  ^-  ,  where  p  is  an  integer. '^  -^  iiut  if  p  is  odd,  tnen  it  is 

P 
easily  seen  tnttt  S(r)  v/ill  have  at  least  two  points  in  li.      To  establish  this 

we  need  only  consider  successive  reflections  in  the  two  sides  which  make  the 
angle  ^  .   Therefore  only  those  fundamental  domains  of  discrete  groups  of  re- 
flections for  which  the  angle  between  any  two  bounding  plajies  is  of  the  form 

~  ,  with  p  integral,  can  be  admissible  regions  for  the  image  method.   Since 

^  [11 

all  the  fundamental  domains  of  discrete  groups  of  reflections  are  known, 

it  cfJi  be  immediately  verified  that  in  two  and  three  dimensions  all  those 

domains  vhich  satisfy  this  additional  condition  are  adm.issible  regions  for  the 

image  method. 

Thus  all  the  admissible  regions  for  the  image  method  have  been  determined. 

Those  in  two  and  three  dimensions  are  shown  in  figures  1  and  2  respectively, 

Somraerfeld  [2]  mentions  most  of  these  ahmissible  regions  in  his  discussion  of 

the  image  method, 

I II ,  Other  boundary  Conditions 

If  the  condition  that  the  normal  derivative  must  vanish  is  replaced  by 
the  condition  that  the  fvmction  vanish  on  some  or  all  of  the  boiandine  Tlanes 
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Figure  1 

Admissible  regions  for  application  of 
image  method;  two  dimensions 


ONE    BOUNDARY 


HALF    SPACE 


-  u  - 

TWO       BOUNDARIES 


WEDGE  (p) 


SLAB 


DUN  D A  R  I  E  S 


\ 


L 


HALF  SLAB        INFINITE      TRIANGULAR     TRIHEDRAL  ANGLE 

PRISMS  (  2.2, p), (2,3,3 )T(2.3,4)* 

(2,3,5)* 

FOUR         BOUNDARIES 


::;> 


INFINITE  RECT.    QUARTER     HALF    TRIANGULAR    PRISMS         PIE   SLICE 


PRISM 


SLAB  (3,3.3)* 

(2,2,2,2) 


(2,3,6)  (2,4,4) 


FIVE    BOUNDARIES 


TETRAHEDRA  HALF    INFINITE  FIN  TE    RIGHT 

(0123)"  (0023)*(0033)*RECTANGULAR  PRISM     TRIANGULAR 

SIX  BOUNDARIES 


RECTANGULAR    PARALLELIPIPED 


Figure  2 


Admissible  regions  for  application  of  image  method; 
three  dimensions. 

For  regions  denoted  by '^^ Laplace's  equation  is  not 
solvable  by  separation  of  variables. 
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of  D,  then  the  image  method  previously  described  mvist  be  modified.   This 
modification  can  be  made  if  it  is  possible  to  assign  a  plus  or  minus  sign 
to  each  point  in  S(r)  in  such  a  way  that  two  points  which  are  images  of  each 
other  in  the  bounding  plane  B  have  the  same  or  opposite  signs  according  as 
the  normal  derivative  or  the  function  vanishes  on  B,  For  then  we  define  G  by 


(2)         G( 


r.r')  =  y     (1)  G*(r".r'). 
?"'cS(r) 


The  sign  of  each  term  in  the  sum  is  that  assigned  to  the  point  r"  in  S(r), 
assuming  that  r  itself  has  a  plus  sign.   It  is  clear  that  (2)  yields  a 
fiinction  which  is  odd  with  respect  to  each  plane  on  which  the  boundary  condi- 
tion reqiiires  that  the  function  vanish,  and  thus  this  condition  is  fulfilled, 
tiimilarly  (2)  is  even  with  rescect  to  every  plane  on  which  it  is  required 
that  the  normal  derivative  vanish,  and  therefore  this  condition  is  also  met. 
Thus  we  see  that  (2)  yields  that  Green's  function  for  the  operator  L  in 
the  region  D  which  satisfies  the  imposed  boiuidary  conditions,  provided  that 
the  previous  conditions  on  L  and  D  are  satisfied,  and  provided  also  that  it 
is  possible  to  assign  plus  and  minus  signs  to  S(r)  in  the  requisite  manner. 
The  conditions  on  L  have  not  yet  been  examined,  but  all  the  admissible  regions 
D  have  been  determined.  We  now  consider  the  possibility  of  assigning  plus 
and  minus  signs. 

Each  region  D  is  a  fundamental  domain  of  a  discrete  group  of  reflections 
which  satisfies  the  extra  condition  that  the  angle  between  any  two  planes 
w^  and  w^  is  of  the  form  -   ,   Now  It  is  a  known  theorem  that  the  discrete  group 

is  completely  characterized  by  the  relations 


(3)        R?   =  1,   (R^  fi  /^-^  =  1 


(When  two  walls  are  parallel,  p^   is  infinite  and  the  corresponding  relation 

is  to  be  disregarded,)   We  make  use  of  this  result  by  supposing  that  the 

boundary  condition  G  =  0  is  imposed  on  the  planes  w,  ..,  w^,  and  5-^  =  0  is 

1      o      On 

iDiposed  on  the  remaining  planes  w<^^^  ...  w^^  .   Now  we  consider  the  group 

generated  by  -R^,    ...,  -B^,  flg^^^ ,  ...  H^.   We  assign  a  plus  or  minus  sign  to 

each  image  according  as  it  is  obtained  from  r  by  a  reflection  bearing  a  plus 
or  minus  sign. 
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In  order  that    this   assignment   be  \inique   it    is   obviously  necessary  and 
siifficlent,   in  virtue  of   the  above    theorem   (see   eq,    (3))»    that 

(4)  (fi^  fij)    ^J  =  1. 

Here  R'  and  R'  stand  for  plus  or  minus  R  and  R  respectively,  according  as 

R.  and  R.  appear  with  plus  or  minus  signs  in  the  generators  of  the  group. 

Thus  we  arrive  at  the  result  that  the  desired  Green's  function  can  be  con- 
structed by  tne  image  method  if  and  only  if  p   is  even  for  every  pair  of 
walls  on  which  different  boundary  conditions  are  imposed.  The  combinations 
of  boundary  conditions  which  can  be  employed  with  each  admissible  region  are 
thus  completely  determined.   In  particular,  if  the  function  is  to  vanish  on 
every  bounding  plane  then  no  additional  condition  on  the  angles  is  required. 

It  is  also  possible  to  extend  the  image  method  to  cases  in  which  the 
boundary  condition  on  some  or  all  bounding  planes  is 

(5)  i  -  »<>  =  0. 

Here  h  is  a  constant  for  each  plane  and  may  be  different  for  different  plajies. 
When  (5)  holds,  we  must  introduce  line,  surface,  and  volume  images,  and  so  on 
for  higher  dimensions.  First  we  will  give  the  result  for  the  case  of  a  single 
plane,  which  we  take  to  be  x '  =0,   Omitting  all  variables  except  x,  we  have 

(6)  G(i,x»)  =  G*(x,x')  +  G*(-x,x«)  +  2h  r   e^^''^'?^G*(T,x')d*| 

It  is  easily  verified  that  (6)  satisfies  (5)  on  x'  =0  provided  that 
aG*(r,x')     3G*  (x.x»). 

(7)         3?  '  -  ax' 

Thus    (6)   is   the  Green's  function  for  the  half -space  bounded  by  the  plane 
x'    =  0  and  satisfies    (5)   on  this   plane. 

To  obtain  tne   Green's  function  for  any  region  D  bovmded  by  planes   on 

aG      Sf 

each  of  which  G  =  0  or  ■r-  =  0  or  t—  +  hG  =  0,  we  proceed  as  follows.  To 

3G  _  « 
each  bounding  plane  we  assign  a  generalized  reflection  operator.  If  c—  -  0 

on  a  plane,  the  corresponding  operator  is  the  ordinary  reflection  operator. 
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while  if  G  =  0  on  a  plane,  the  operator  is  the  ordinary  reflection  operator 
multiplied  by  -1.   However,  if  the  condition  (5)  holds  on  a  plane,  the  corres- 
ponding operator  assigns  to  each  point  its  ordinary  image  pl\i6  a  line  of  images 
on  the  half -ray  vhich  is  normal  to  the  plane  and  extends  from  the  image  point 
away  from  the  plane.  The  ordinary  image  is  assigned  the  same  strength  as  the 
antecedent  point,   Jiiach  point  on  the  ray  is  assigned  a  strength  equal  to  2he 
times  the  strength  of  the  antecedent  point.  The  strength  of  the  original 
point  r  is  +1,  tiutt  of  an  image  in  a  plane  on  which  ~-  =  0  is  also  +1,  and  that 
in  a  plane  on  which  G  =  0  is  -1,   In  this  way  all  the  image  points  and  their 
strengths  are  defined,   Now  the  Green's  function  G(r,r*)  for  the  region  D  is 

constructed  by  forming  the  sum  J]  G*(x",x')  A(r"),  where  the  sum  extends  over  all 

r" 
images  of  r  and  where  A(r" )  is  the  strenpth  of  each  imaf e ,  The  sum  will,  of  course 

include  line,  surface,  and  volume  inteprals,  etc. 

In  order  for  this  construction  to  succeed  it  is  again  necessary  that  no 

image  lie  within  L,  except  r  itself.   Therefore  I)   must  be  one  of  the  regions 

previously  determined.  If  the  boundary  condition  on  some  planes  is  that  G  =  0, 

then  the  angles  between  these  planes  and  the  others  must  be  even  submultiples 

of  TT.   The  'other'  planes  are  those  on  wnich  ~~     =  0  or  ~-  +  hG  =0,  Further- 

Bn        dn 

more,  in  order  that  no  line  of  images  lie  in  D  it  is  necessary  that  the  angle 
between  any  plane  on  which  (5)  is  imposed  and  any  other  plane  be  ~  ,  as  is 
easily  verified. 

These  results  completely  characterize  the  regions  and  combinations  of 
boundary  conditions  to  which  the  extended  image  method  applies, 

I V .   Mmissible  Differential  Equations 

In  order  that  the  image  method  apply,  the  only  reo^uirement  on  the  linear 
differential  operator  L  is  that  each  term  in  the  sura  (l),  (2),  or  in  the  stun 
described  at  the  end  of  the  previous  section  also  be  a  solution  of  MG  =  0, 
The  necessary  and  sufficient  condition  for  this  to  be  the  case  is  that  M  admit 
the  group  of  reflections  in  the  boundaries  of  D,  i,e,,  that  M  be  invsu-iant 
under  these  reflections.   In  the  case  of  the  mixed  boundary  condition,  (7)  must 
also  hold.   These  considerations  determine  the  differential  equations  to  which 
the  image  method  applies. 


t  Iteration  of  the  generalized  reflection  operation  is  defined  to  be  the 
identity. 
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As  an  eiaarple,   let   us   consider   operators   M  in  two  variables  with 
constant  coefficients.      The   eqtiation  Mu  =  0   is   equivalent   to 

(S)       au^  +  bu^y  +  cu^  +  du^  +  eu^  +  fu  =  0. 

If  x',  y'  is  the  image  of  x,y  in  a  line  making  the  angle  »  with  the  x-axis, 
tnen  (8)  becomes  in  the  x',  y'  coordinates 

(9)  (a  co3^2»  +  b  sin  26  cos  26  +  c  8in^26)  \,j^, 
+  (a  sin^  26  -  b  sin  26  cos  26  +  c  co8^26)u  ,  , 
+  (2  a  Bin  26  cos  26  +  b  sin^  26  -  b  co8^26  -  2c  sin  26  cos  26)u  ,  , 

*  y 

+   (d  cos  26  +  e  sin  26)u  ,     +   (d  sin  26  -  e  cos  26)u   ,  +  f u  =  0 . 

In  order  that  (8)  be  invariant  under  this  reflection,  (9)  must  coincide  with 
(8),  Therefore  the  coefficients  in  (9)  must  be  multiples  of  the  correspond- 
ing coefficients   in   (8),      If  f   is   not    zero,    then  the  multiple   is   one. 

We   th'os    obtain  e.  set   of  five   linear   homof^sneous   equations   for  the   coef- 
ficients a  to  e,   which  separate  into  a  set  of  three   for  a,b,c  aiid  a  set  of 
two  for  d,   e.      The   ranic  of   eacli  coefficient   matrix   in  these  eq'jations    is   one, 
so   that   the  first   three  equations   have   a  two-parameter  set   of   solutions,   and 
the   last   two  a  one-parameter  set.      They  are   given  by 

(10)  -7-  =  tan  26 
a-c 

(11)  -£    =       si  a  ^» 

d         1   -  cos   26 

If  both  b  =  0  and  a-c  =0,  (lO)  is  unnecessary  and  the  first  three  eqioations 
for  a,  b,  c  hold  for  every  value  of  6, 

The  first  result  is  then,  that  if  f  i  0,    {8)    is  invariant  under  reflection 
only  for  those  angles  6  which  satisfy  (10)  and  (ll).  If  both  b  and  a-c  are 
zero,  (10)  Is  unnecessary,  while  if  both  d  and  e  are  zero,  (ll)  is  \annecessary. 

In  the  interval  0  ^  6  <  tt,  (ll)  has  only  one  solution  and  (10)  has  only 
two  which  differ  by  -  ,   Therefore  when  (10),  (ll)  apply  there  is  at  most  one 
angle  6  satisfying  Doth.   If  (ll)  doesn't  apply  but  (lO)  does,  there  are 
exactly  two  angles,  differing  by  ^,  which  satisfy  (lO),   If  neither  (lO)  nor 
(11)  applies  then  6  is  unrestricted.   Thus  we  have  the  following  results  pro- 
vided f  t  0. 


e  =  cL  =  b=0,a  =  c  all© 

e,  d  not  botn  zero,  b  ^  0,  a  =  c  one  0  given  by  (ll) 

e,  d  not  both  zero;  b,  a  -  c  not  both  zero   at  most  one  ^   given  by  (10,  ll) 

e=d=0;  b,  a-c  not  ooth  zero  two  angles  ft,  »  +  -  given  by  (lO)  , 

If  f  =0,  then  in  addition  to  the  above  results  we  find  that  if  a  =  -  c  other 
possible  values  of  ©  are  roots  of 

2a    sin  2©  =  i 


tan  2©   1.   •   1  ^   „  OQ      ^ 
b     1  +  cos  2©      e 


TT 


iigain  the  first  has  two  roots  differing  by  "  and  the  second  only  one  root. 
Tqus  if  both  apply  there  is  only  one  value  of  ©,  and  if  only  the  first  applies 
there  are  two  values  of  ©  differing  by  "  . 

V,  Conclusion 

The  image  method  has  been  defined  for  regions  in  ihjclidean  n-space 
bounded  by  planes,  and  all  the  regions  to  which  it  applies  have  been  deter- 
mined for  every  choice  of  boundary  conditions.   The  bo'indary  conditions  con- 

sidered  are  G  =  0,  ^  =0,  and  r~  +  hG  =  0.  For  each  region  the  admissible 

on         On 

differential  equations  have  also  been  characterized.   It  is  interesting  to 
note  that  of  the  problems  which  can  be  treated  by  the  image  method  not  all  can 
be  solved  by  separation  of  variables. 

A  generalization  of  the  image  method,  namely  inversion,  applies  to  regions 
bounded  by  spheres  but  only  worxs  for  Laplace's  equation.   It  would  be  inter- 
esting to  determine  all  the  regions  to  which  this  method  applies. 

The  image  method  can  also  be  applied  to  differential  equations  on  the 
surface  of  a  sphere  and  to  regions  bounded  by  great  circles.   The  present 
results  yield  all  the  regions  on  the  sphere  to  which  this  method  applies.  We 
need  only  consider  tnose  Euclidean  regions  with  one  vertex,  and  for  each  con- 
struct a  sphere  with  center  at  this  vertex.  The  region  cut  out  on  the  spheri- 
cal surface  by  the  Euclidean  region  is  an  admissible  image  region  on  the  sur- 
face of  the  sphere,  and  conversely. 

The  method  can  also  be  extended  to  other  surfaces,  such  as  Riemann  sur- 
faces, as  was  done  by  Sommerfeld.   In  this  case,  however,  the  main  difficulty 
is  the  constrtiction  of  the  Green's  ftmction  for  the  full  space. 

Another  extension  is  to  regions  bounded  by  an  even  number  of  pairwise 
parallel  planes  (or  other  surfaces)  with  the  boundary  condition  of  periodicity. 
This  conaition  is  that  the  boundary  values  of  the  solution  miust  be  the  same 
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at  corresponding  points  on  two  parallel  congruent  boundary  planes,  'Jn.e 
admissible  regions  can  be  found  by  considering  discrete  groups  of  transla- 
tions.  If  periodic  boundary  conditions  are  inroosed  on  some  ul^nes .  and  G-  =  0 
an  ~      others,  then  more  general  discrete  groups,  involving  trans la- 


or 


tions   and  reflections,   will  play  a  role.      If  periodic  boundary  conlitions   are 
imposed  on  two  planes   which  meet   at  an  angle,    then  discrete  groups    involving 
rotations   must   be   considered, 

i"or  certain  regions   only  a  finite  number  of    images    occur,   and   ther'=fore 
the   sums   used   in  defining  the  Green's   function  are  finite.      In  otner  cases 
the  convergence   of   the  sum  must   be   examined,      Tnis   has   been  done   in  some 
cases.    [2] 

It  is  a  pleasure   to  thank  Professor  VJilhelni  Magnus  Jjr  his  helpful 
advice   and  criticism  in  connection  with  this  paper. 
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